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Abstract. Two-phase flow calculations are presented to investigate the thermo-hydraulical effects of the 
interaction between 300 kj proton pulses (2 ms long, 1.3 GeV) with a closed mercury loop which can 
be considered as a model system of the target of the planned European Spallation Source(ESS) facility. 
The two-fluid model consists of six first-order partial differential equations that present one dimensional 
mass momentum and energy balances for mercury vapour and liquid phases are capable to describe quick 
transients like cavitation effects or shock waves. The absorption of the proton beam is represented as 
instantaneous heat source in the energy balance equations. Densities and internal energies of the mercury 
liquid-vapour system are calculated from van der Waals equation, and the general method how to obtain 
such properties is also presented. A second order accurate high-resolution shock-capturing numerical scheme 
is applied with different kinds of limiters in the numerical calculations. Our analysis show that even 75 
degree temperature heat shocks cannot cause considerable cavitation effects in mercury in 2ms long pulses 



PACS. 47.55Bx Cavitation, - 47.55Kf Multiphase and particle-laden flows - 
dynamics 



47.90a-f Other topics in fluid 



1 Introduction 

One non-destructive material research method is neutron 
scattering. Free neutrons for neutron beams for research 
purposes need to be extracted from their bound states of 
atomic nuclei. Energetic neutron beams can be produced 
in fission of heavy elements (e.g. "^^^U) or by spallation. 
In fission of ^35^ ^90 MeV heat is released for each ex- 
tracted fast neutron while in the spallation process only 
about 30 MeV heat is deposited per fast neutron. The 
deposited heat has to be removed by cooling and it ulti- 
mately becomes a limiting thermodynamic factor for the 
amount of neutrons produced. As a second distinct advan- 
tage of pulsed spallation sources over continuous ones is 
that a larger part of the neutrons produced can be deliv- 
ered to the sample in monochromatic beams. These two 
advantages of spallation sources make is possible to con- 
struct more powerful neutron sources with larger neutron 
flux than ever before. The simple goal of the planed Eu- 
ropean Spallation Source(ESS) is to provide Europe with 
the most powerful neutron facility. A choice of a 5 MW 
proton beam power at 1.3 GeV proton energy with 111 mA 
proton beam current and with 16.66 Hz repetition rate of 

2 ms long neutron pulses will produce time average ther- 



mal neutron flux density of 3.1 x / cm? s in the ESS 

mercury target. The proton pulse causes a thermal and a 
pressure shock in the target which may cause cavitation 
or tensile stress. The question of cavitation erosion [T] has 
crucial importance in the constructional planning of any 
spallation neutron source target facility. A detailed anal- 
ysis of the planned ESS can be found elsewhere [2]. In the 
following study we present and analyze a one-dimensional 
six-equation two-fluid model which is capable to describe 
transients like pressure waves, quick vapour void fraction 
creation and annihilation which is proportional to cavi- 
tation caused by energetic proton interaction in mercury 
target. 

Our model has a delicate numerical scheme and capa- 
ble to capture shock waves and describe transient waves 
which may propagate quicker than the local speed of sound 
3J . Most of the two-phase models have numerical methods 
which describes usual flow velocities. 

Our model can successfully reproduce the experimen- 
tal data of different one- or two-phase flow problems such 
as ideal gas Riemann problem, critical flow of ideal gas in 
convergent-divergent nozzle, column separation or cavita- 
tion induced water hammer, rapid depressurization of hot 
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liquid from horizontal pipes or even steam condensation 
induced water hammer 4]. 

According to our knowledge there is no real two-phase 
flow calculation for mercury flow system. Some timorous 
attempts were presented with the help of some commercial 
three dimensional industrial codes like Fluent or ANSYS 
[5ll6] but the results are questionable. Some results show 
complete and immediate vaporization during the first pro- 
ton pulse, which is contradictory to experimental observa- 
tions. There are some studies for three dimensional numer- 
ical simulation of magnetohydrodynamic processes in the 
muon colliders mercury target. These studies take strong 
external magnetic fields into account [T but consider sin- 
gle phase only, neglecting vaporization or condensations. 
The liquid phase of mercury was modeled using the stiff- 
ened politropic equation of state and the vapour phase was 
considered to be ideal gas. There is a literature survey on 
various fluid flow data for mercury from the politropic 
equation state [8] which can be directly applied in calcu- 
lations. There are also different equations of state (EOS) 
available for mercury from microscopic molecular simula- 
tion [QlfTO] of from macroscopic theories like virial expan- 
sion [TT] or from generalized van der Waals equation like 
the Redlich-Kwong equation pjj or the like [13^ . Thermo- 
dynamical and flow properties of other liquid metals are 
also in the focus of recent scientific interest p^lT2] . 

In the next sections we introduce our applied model, 
give a detailed analysis about phase transitions and present 
pressure wave results. 



2 Theory 
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Index 1 refers to the liquid phase and index g to the 
gas phase. Nomenclature and variables are described in 
Table I. Left hand side of the equations contain the terms 
with temporal and spatial derivatives. Hyperbolicity of 
the equation system is ensured with the virtual mass term 
CVM and with the interfacial term (terms with pi). Terms 
on the right hand side are terms describing the inter-phase 
heat, mass (terms with Fg vapour generation rate) volu- 
metric heat fluxes Qig^ momentum transfer (terms with 
Ci), wall friction Fg^aii, and gravity terms. A detailed 
analysis of the source terms can be found in [3]. The last 
term in the energy equations Eipuise{x^t) represents the 
deposited energy from the proton beam and will be spec- 
ified later on. 

Two additional equations of state (EOS) are needed to 
close the system of equations (Eq. 1-6) 
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2.1 Theory of two-phase flow 

There is a large number of different two-phase flow mod- 
els with different levels of complexity [i5,16_ which are all 
based on gas dynamics and shock-wave theory. In the fol- 
lowing we present our one dimensional six-equation equal- 
pressure two-fluid model. The density, momentum and en- 
ergy balance equations for both phases are the following: 
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Partial derivatives in Eq. 7 are expressed using pressure 
and specific internal energy as an input. In the following 
we show how the liquid-steam table - a sixfold numerical 
table - (p, T, pi,ui, pg,Ug) can be created for mercury from 
an arbitrary EOS. To avoid technical difficulties we do 
not modify (Eq. 1-6) including the used analytic EOS, 
just create a numerical liquid-steam table. In this manner 
arbitrary two phase-flow systems can be investigated with 
the same model in the future (e.g. lead-bismuth eutectic, 
liquid Li or He). Liquid metal systems can operate on 
low (some bar) pressure and have larger heat conductivity 
than water which can radically enhance thermal efficiency. 

We start with the usually parameterized van der Waals 
EOS from 

RT a 

where R=8.314 J/mol/K is the universal gas constant and 
parameters a and b are related with the critical molar vol- 
ume (Vc) temperature (Tc) and pressure {pc) of the con- 
sidered fluid: a = 9PcV^ b = Vc/3 For the critical temper- 
ature and pressure of Hg the Tq = 1733 ± 50K and pc = 
160.8 ± 5MPa data were taken from [T^]. T, p and V are 
the temperature the pressure and the volume, respectively. 
(We mention than in [12| the parameter a(= QPcV^) is un- 
corectly given.) The fluid density with the corresponding 
saturated vapour density can be easily determined from 
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the EOS with the well-known Maxwell construction. To 
obtain the internal energies for both phases is a bit more 
difficult task. We start with the second law of thermody- 
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where s is the entropy and u is the internal energy. With 
the Maxwell relations (f^)g = — (w)y with 
the following working equation 
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The internal energy is a thermodynamical potential there- 
fore the choice of the zero point can be defined arbitrary, 
we took T = 253.1AK which is 10 degree above the melt- 
ing point of solid mercury at normal pressure. The heat 
capacity at constant volume Cy may in turn be calculated 
from the heat capacity at constant pressure Cpwith the 
thermodynamic relation 



(11) 



where = Q^t is the thermal expansion coefficients. 

(To avoid further misunderstanding in this study we use 
ax for the thermal expansion coefficient and a for vapour 
void fraction.) Polynomial fits for the temperature depen- 
dence of experimental data of heat capacity Cp and ex- 
pansion coefficient ar [12](or 0) help us to calculate the 
internal energy of the liquid state. Finally, the internal en- 
ergy of the corresponding gas phase has to be determined. 
The critical temperate of mercury is at Tc = 1733 ± 50K. 
In the temperature range of 270-500 K (which is our recent 
interest) the experimental heats of vapourization data [10| 
can be satisfactory fitted with a linear function. With this 
method a two-phase liquid-steam table was constructed 
between 270-500 Kelvin of temperatures and 7 to 7 • 10^ 
Pascal pressure. 

Additional flow properties of mercury like dynamic 
viscosity and heat transfer coefficients are approximated 
with piecewise continuous temperature dependent func- 
tions from f8]. The surface tension was considered as a 
linear function of temperature [18;. 

The effect of the 300 kJ proton pulse was treated as a 
sudden thermal shock which means an additional source 
terms in both energy equations Ei_puise{x, t). The deposited 
energy of the proton beam in the mercury target is propor- 
tional to the density. With the introduction of the mixture 
density pm — apg + (1 — '^)Pf the interaction between the 
proton-mercury two-phase fiow can be further improved. 
According to experimental proton beam analysis the spa- 
tial energy distribution perpendicular to the propagation 
direction has a parabolic shape [l9l[20] with a diameter 
of 20 cm. To describe well-defined finite duration we use 
sin^ envelope with r = 2ms. 
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The effective range of 1.3 GeV protons in mercury can be 
calculated with the Bragg theory and gives 41 cm |21j . 
Experimental consideration state that 47 percent of the 
original 5 MW beam power is absorbed in the target which 
is 2.37 MW. In the planned ESS facility a train of 16.66 
proton pulses will come with 2ms long pulse duration and 
the total sum of these pulses give the 5 MW beam power. 
Hence, the peak power parameter Eq has to be normal- 
ized in such a way that the spatial and time integral of 
Ei^puise{x,t) gives the absorbed 2.37 MW power of the 
original 5 MW beam. The system of Eqs. (1-6) represents 
the conservation laws and can be formulated in the fol- 
lowing vectorial form 



A— B— S 

dt dx 



(14) 



where tf^ represents the vector of the independent non- 
conservative variables ^{p^a,Vf,Vg,Uf,Ug), A and B are 
the matrices of the system and S is the source vector of 
non-differential terms. These three quantities A, B and S 
can be obtained from Eq. (1-6) with some algebraic ma- 
nipulation. In this case the system eigenvalues which rep- 
resent wave propagation velocities are given by the deter- 
minant (iei(B — AA). An improved characteristic upwind 
discretization method is used to solve the hyperbolic equa- 
tion system (|14p . The problem is solved with the combina- 
tion of the first- and second-order accurate discretization 
scheme by the so-called flux limiters to avoid numerical 
dissipation and unwanted oscillations which appear in the 
vicinity of the non-smooth solutions. Exhaustive details 
about the numerical scheme can be found in 4,3'. 



2.2 Liquid-vapour phase transition in the metastable 
region 

Water boils at 100 Celsius (373.15 K) under atmospheric 
pressure; this is a well-known, but not entirely correct 
piece of the common knowledge. Boiling is usually de- 
fined (at least phenomenologically) when liquid-vapour 
phase transition happens not only at the already exist- 
ing interfaces, but within the bulk liquid too. For water, 
it happens usually at the already mentioned 100 Celsius, 
but not always. Overheating of liquids is a phenomenon, 
known for every chemistry students; one can exceed the 
boiling point with a few degrees, without getting boil- 
ing, but then it can happen suddenly, exploding the whole 
amount of liquid (and often the container too) [Hj . In the 
following, we are going to explain this phenomenon and 
show its importance in the cavitation of mercury. Liquid 
can co-exists with the vapour of the same material, with- 
out any problem. The conditions (temperature and pres- 
sure) where they co-exist are described the vapour pres- 
sure curve (also called saturation or co-existence curve). 
Liquid and vapour states can be described by EOS; like 
van der Waals. A schematic isotherm (describing pressure 
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and volume on a constant temperature) can be seen on 
Figure 1/a. The isotherm has two extrema (marked as 
B and D), these are the so-called spinodal points (liquid- 
vapour and vapour- liquid; LV and VL spinodals). Between 
the two spinodals, the system would be unstable, due to 
the negative compressibility, therefore these states (the 
ones on the curve between points B and D) cannot exist. 
The equilibrium conditions can be calculated by using the 
Maxwell construction: a line (parallel to the V-axis) has 
to be drawn in a way that the area between the isotherm 
and the Maxwell-line between points A and C and C and 
E has to be equal. Then the intersects (A and E) gives the 
co-existing liquid and vapour volumes (or densities) and 
the equilibrium pressure on the given temperature. Plot- 
ting the pressures on different temperatures, one would 
obtain the vapour pressure curve, like the solid line on 
Figure 1/b. It can be seen that points A and E are not 
special points of the isotherm. The liquid is not forced to 
boil at point A; it would be forced only at point B (where 
liquid phase cannot exists any more). Plotting Bs at dif- 
ferent temperature, one would obtain the so-called liquid- 
vapour spinodal (dashed line on figure 1/b), the stabil- 
ity limit of liquid state. The vapour-liquid spinodal (dot- 
dash line on figure 1/b) is important when we have over- 
saturated vapour; we are going to neglect it here. The AB 
and DE parts on the isotherm are metastable; on Figure 
1 /b these parts are represented by the region between the 
vapour pressure curve and the LV spinodal. Liquid (with- 
out co-existing vapour phase) can exist in this region; even 
can exist under negative pressure P^[f25] The real boiling 
happens in this metastable region. Close to the vapour 
pressure curve the liquid is only slightly metastable, can 
live for long time without nucleating vapour bubbles; far 
away from it (close to the spinodal) the liquid will be 
very metastable and cavitate (boil) with higher proba- 
bility. The bubble nucleation can happen in two different 
ways. The heterogeneous nucleation happens when the liq- 
uid already has a nucleus, usually a tiny bubble hidden in 
a crevice of the wall or stuck onto a floating particle. Due 
to the small size (i.e. high curvature) the micro-bubble can 
be in equilibrium with a metastable liquid for a while, but 
when the temperature is too high or the pressure is too 
low, it will initiate boiling. The other process is the homo- 
geneous nucleation. In that case, the initial micro-bubble 
will be produced by the density fluctuations of the liquid; 
when the fluctuation is big enough to call it "bubble", 
then it will initiate the boiling. In everyday life, boiling 
happens by heterogeneous nucleation, practically in the 
immediate vicinity of the vapour pressure curve. In clean 
liquids (like distilled water) the boiling can happen much 
farther. It is a well-known practice to avoid overheating 
(and the explosion-like boiling, following it) to put some 
bubble seed into the liquid, like a few pieces of sponge- 
like pumice (or boiling-) stones. In these nucleation pro- 
cesses -especially in homogeneous nucleation - time is also 
an important factor; a liquid can endure high overheat- 
ing/stretching for a small period of time [22 j l23 ll24i 25 . 
Therefore one cannot draw a well-defined line as nucle- 
ation limit onto the phase diagram (Figure 1); the line 



drawn by us is only for demonstration, the exact location 
depends on the purity of liquid, the amount of external 
disturbances (even cosmic rays can generate bubbles in 
metastable liquids) and - in a great extent - on the time 
scale. On Figure 1/c, a magnified part of Figure 1/b (with- 
out the VL spinodal, which is irrelevant in our case) can be 
seen. K marks a state, where the liquid is in stable liquid 
phase; there is no vapour phase present. To obtain phase 
transition, the temperature can be increased or the pres- 
sure has to be decreased. By increasing the temperature 
(and keeping a constant pressure), the vapour pressure 
curve will be reached at point L. This is the first point, 
where the liquid can boil and vapour phase might appear, 
but in clear and undisturbed liquid, the probability of boil- 
ing here is very small. Increasing the temperature further, 
the nucleation limit will be reached (point M); here the 
phase transition will surely happen, due to heterogeneous 
or homogeneous nucleation, forming initially small, but 
continuously growing separated bubbles, which later can 
merge into a continuous vapour phase. When the liquid 
is perfectly clean, all disturbances are suppressed and the 
heating is very fast, etc., then this nucleation limit can 
be pushed very close to spinodal limit (point N), where 
liquid phase cannot exist any more. When the phase tran- 
sition happens at the spinodal limit, one will obtain two 
bi-continuous phases (liquid and vapour) , instead of a con- 
tinuous (liquid) and an separated (vapour bubbles) one, 
obtained during nucleation. We should remark, that for 
the first appearance of the vapour phase, the system will 
jump back to the vapour pressure curve (which will be 
detected as a pressure jump). Changing the pressure at 
constant temperature, one would reach the vapour pres- 
sure curve at point O, then the nucleation limit at point 
P, finally the spinodal limit at point Q, with the same 
results as it happens with temperature increase. To see 
the extent of the effect, we will give numerical examples 
for water and for mercury. For water, starting from room 
temperature (293.15 K, 1 bar), we will reach the vapour 
pressure curve 373.15 K. Increasing the temperature, boil- 
ing might happen any time; the highest experimentally 
obtained value for overheating (i.e. point M) was around 
570 K [27] giving almost 200 K overheat. The spinodal 
temperature for water on atmospheric pressure is still de- 
bated, it has to be located above the previously mentioned 
overheating limit, but certainly below the critical temper- 
ature of 1733 ± 50 K. Also for water, by decreasing the 
pressure, the vapour pressure curve (point O) would be 
reached at 0.025 bar pressure. The experimental limit of 
stretching is -1200 bar [27], where the estimated spinodal 
(depending on the model) is between -2000 and -4000 bar. 
For mercury at 7 bar (which is the working pressure for 
the mercury in the ESS) the boiling point is at 760 K, it 
is very far from the working temperature (which is close 
to room temperature, 300 K). The limit of overheating 
is not known, but surely below the critical temperature, 
which is around 1700 K. Concerning pressure drop, the 
vapour pressure of mercury around room temperature is 
almost zero (less than 2*10-6 bar); concerning the fact 
that the working pressure is 7 bar, the possibility of a 
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Fig. 2. The schematic geometrical model of the ESS target 
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Fig. 1. A subcritical isotherm of a van der Waals type fluid 



pressure drop of this extent is very improbable. The mea- 
sured nucleation limit of mercury at room temperature 
is in the -2 to -425 bar range depending on purity[26,; 
therefore to get bubbles, the pressure should drop from 
+7 bar to [-2 bar for a longer period. The absolute (spin- 
odal) limit is unknown. Although in the ESS, pressure de- 
crease and temperature increase happens simultaneously, 
the working conditions are so far from the beginning of 
the boiling region (vapour pressure curve) that the pos- 
sibility to reach it is negligible, except under special cir- 
cumstances. First, there is a possibility for fast pressure 
oscillation after the proton pulse; the amplitude can be 
even 300 bar [29] , which would be enough to cause cavita- 
tion. The other scenario would require gas-contamination 
(pre-existing gas bubbles in the mercury); in that case 
even a tiny pressure decrease or temperature increase can 
cause the growing of these micro-bubbles, mimicking boil- 
ing [29j . Non-uniform temperature and pressure distribu- 
tion can cause shear stresses, which can also cause cav- 
itation in the liquid. Finally, the proton beam itself can 
initiate cavitation, but only when the metastable states 
are already reached. We can conclude, that with a few 
bar pressure drop and a few tens of K temperature in- 
crease, the cavitation in the pure mercury target has low 
possibility. On the other hand, concerning the reported 
cases of cavitation in similar facilities 28J indicate, that 
either the conditions (T,p) might change more drastically 
or some phenomena, neglected by us (like pressure oscil- 
lation, shear stresses, etc.) can play more important role. 



3 Results and discussion 

The ESS mercury target loop is a complex facility with 
various pumps, heat exchangers and tanks [5]. We model 
however with a simple six-sided closed loop (see Figure 



1.) of a pipe with diameter of 5 cm and total length of 5 
m. The original temperature of the mercury is T = 300 
K with pressure of 7 bar and flow velocity of t; = 4.6 
m/s. The proton beam interacts with a mercury via a 
20 X 5 cm^ window. A simple calculation shows that (47 
% of 300 kJ =) 141 kJ of energy will heat up 10 kg of 
mercury. The temperature jump of At — 75K is expected 
for ESS proton beam pulses. We applied a single pulse 
shot at time equal to zero and propagated Eqs. (fT4|) to 
tmax = 4 * 10~^ sec. A second order numerical scheme 
was used with the MINMOD flux limiter [3]. For a satis- 
factory convergence the Courant number which measures 
the relative wave propagation speeds of the exact solution 
and the numerical solutions was set to CFL — 0.6. The 
pressure history of the beam-target interaction point is 
presented in Fig 2. After the initial pulse at t = 1.6ms 
the pressure reaches its maximal value which is 50 per- 
cent higher than the original pressure. After the pulse the 
pressure does not fall below the initial pressure and the 
temperature will cool down to 300 K. The mercury vapour 
void fraction was originally set to zero (a = 10"^^) which 
did not chang during the time propagation allowing only 
" nanobubbles" , too small to act as cavitation nuclei. If we 
consider one or two percent initial vapor void fraction (as 
a model for small bubbles) than a quick condensation can 
be observed. If we apply an elastic pipe with an elastic- 
ity of 2 X 10^^ N/m"^ Young's modulus (which are usual 
for steel) or / and include or exclude any kind of additional 
wall friction 4 for the fluid the pressure peaks will not 
be changed. This is clear fingerprint that the tube is still 
rigid enough. There is a strong indication that mercury is 
a non-wetting fluid on steel surface so the wall friction is 
negligible. 



4 Summary and outlook 

With the help of a one dimensional two-phase flow model 
we calculated the induced pressure waves and vapour void 
fractions in mercury induced by energetic proton beams. 
Our analysis showed that no vapour bubbles or cavitation 
effects can be seen after the first absorbed proton pulse. 
Further, in depth analysis is in progress to investigate ge- 
ometrical effects of the mercury target loop which is a 
complex facility with various pumps, heat exchangers and 
tanks ..2\. Our model can include abrupt area changes, or 
convergent-divergent pipe cross section changes, or even 
heat exchangers. We modeled however with a simple six- 
sided closed loop (see Figure 1.) of a pipe with diameter 
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Fig. 3. Time liistory of pressure at tlie point of the proton 
impact for the model in Fig. 2. with rigid tube walls. 

of 5 cm and total length of 5 m. The original temper- 
ature of the mercury is T — 300 K with pressure of 7 
bar and flow velocity of w = 4.6 m/s. The proton beam 
interacts with a mercury via a 20 x 5 crri^ window. A sim- 
ple calculation shows that 141k J of energy will heat up 
10 kg of mercury with a, At = 75K. We applied a sin- 
gle pulse shot at time equal to zero and propagated Eqs. 
([HI) to tmax = 4 * 10~^ sec. A second order numerical 
scheme was used with the MINMOD flux limiter ^3J. The 
Courant number which measures the relative wave prop- 
agation speeds of the exact solution and the numerical 
solutions was set to CFL = 0.6. The pressure history of 
the beam-target interaction point is presented in Fig 3. Af- 
ter the initial pulse at i = 1.6ms the pressure reaches its 
maximal value which is 50 percent higher than the original 
pressure. After the pulse the pressure does not fall below 
the initial pressure and the temperature will cool down to 
300 K. The mercury vapour void fraction was originally 
set to zero (a = 10~^^) which did not changed during the 
time propagation allowing only " nanobubbles" , too small 
to act as nucleus for cavitation. The question of the vapor 
void fraction, pipe elasticity or the liquid wall friction was 
examined also. 

We would like to emphasize that further in-depth anal- 
ysis is needed to clear up the question of a long pulse train. 
The question of different equations of state will be inves- 
tigated also. As a long term interest we also planned to 
investigate other liquid metal (e.g. bismuth-lead eutectic 
or liquid lithium) or liquid helium systems which can be 
interesting as a cooling media for new type of nuclear reac- 
tors. Liquid metal systems can operate on low (some bar) 
pressure and have much larger heat conductivity than wa- 
ter which can radically enhance thermal efhciency. 
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Table 1. Nomenclature used in the two-phase flow equa- 

tions(Eq. 1-6. ) 

A pipe cross section (m"^) 

Ci internal friction coefficient {kg/rn^) 

CVM virtual mass term [N/m^) 

Ci specific total energy [e = u + /2] (J/kg) 

Ff,waii wall friction per unit volume (N/m^) 

g gravitational accclcration(m/s^) 

hi specific enthalpy [h = u + p/p] (J/kg) 

p pressure (Pa) 

Pi interfacial pressure pi = pa{l — a) (Pa) 

Qig interf.-liq./gas heat transf. per vol. rate (W/m^) 

t time (s) 

Ui specific internal energy (J/kg) 
Vi velocity (m/s) 

Vr relative velocity [vr = Vg — Vf ]{m/s) 

w pipe velocity in flow direction {m/s) 

X spatial coordinate (m) 

r vapour generation rate (kg/m^) 

a vapour void fraction 

Pi density {kg/m^) 

6 pipe inclination 



